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Metamaterial cloaking has been proposed and studied in recent years following several interesting 
approaches. One of them, the scattering-cancellation technique, or plasmonic cloaking, exploits the 
plasmonic effects of suitably designed thin homogeneous metamaterial covers to drastically suppress 
the scattering of moderately sized objects within specific frequency ranges of interest. Besides its 
&X)' inherent simplicity, this technique also holds the promise of isotropic response and weak polarization 

^ I dependence. Its theory has been applied extensively to symmetrical geometries and canonical 3D 

•^r ■ shapes, but its application to elongated objects has not been explored with the same level of detail. 

We derive here closed-form theoretical formulas for infinite cylinders under arbitrary wave incidence, 
and validate their performance with full-wave numerical simulations, also considering the effects of 
finite lengths and truncation effects in cylindrical objects. In particular, we find that a single 
isotropic (idealized) cloaking layer may successfully suppress the dominant scattering coefficients of 
ri ' moderately thin elongated objects, even for finite lengths comparable with the incident wavelength, 

Oj, providing a weak dependence on the incidence angle. These results may pave the way for application 

I ■ of plasmonic cloaking in a variety of practical scenarios of interest. 

'Tj [ I. INTRODUCTION 

. ^ ■ The a ppli cation of metamaterials to cloaking and invisibility has been explored in several exciting papers in recent 

C/3 ' years |l|-[l5j. In these contributions, several alternative realizations and techniques have been discussed, and the 

^~>| exotic properties of several classes of metamaterials have been shown to be possibly tailored in order to provide much 

reduced scattering from finite-sized objects in different configurations and schemes, for a wide range of frequencies of 

interest. Recent reviews of the various possibilities are available for the interested reader (see e.g. Refs. |16l.ll7t). 

One such possibility takes advantage of the anomalous scattering response of thin plasmonic layers. As shown in 

CN ' Ref. |15| , artificial plasmonic materials with low or negative effective permittivity may provide scattering cancellation 

via their local negative polarizability. This technique, called plasmonic cloaking, is consistent with earlier works that 

have speculated how a composite particle combining positive and negative permittivity may provide identically zero 

^Q ' scattering in the static limit |18| . Recent study of the dynamic case has shown that plasmonic cloaks may suppress 

f»o ■ not only the dominant dipolar scattering from moderately sized objects, but also higher-order multipolar orders for 

larger scatterers [1^. In this vein, it is worth stressing that by "cloaking" we mean strong, or maximized, scattering 

reduction over a finite frequency band - not necessarily complete invisibility, since residual scattering orders may 

p^ ' always make the scatterer detectable to a certain extent. Still, significant reduction of visibility is achievable with the 

proper plasmonic cloak design, as shown in several recent papers [20|-|28[ and as discussed in the following. 

Plasmonic cloaking has been shown to offer several i ntriguin g properties in a variety of setups. Examples include 
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intrinsic robustness to frequency and desi gn variations [20|-|22|, straightforward extension to arbitrary collections of 



objects [23| and multi- frequency operation |2J,|25[. Moreover, the admission of fields inside the cloaked region, peculiar 
Jh [ to this technique, may be used to suppress the inherent scattering from receiving antennas and sensing devices, which 
. . . may open several interesting venues in non- invasive probing and sensing applications [26[ . Extensions to ultrathin 
surface cloaks have also been put forward in this same context [27[. All these studies were conducted, for sirnplicity, 
for a canonical spherical object, illuminated by a generic polarization of the impinging field. Recent studies [28| have 
shown that analogous concepts may be applied to more complex 3D geometries, based on the overall robustness of 
the integral cancellation effect. 
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FIG. 1: Schematic of the cloaked cylindrical object problem and its geometry. A dielectric core of length L and diameter 2a 
is covered uniformly, except the ends, by a metamaterial shell of radial thickness (ttc — a). The object is illuminated by an 
arbitrarily polarized plane wave, incident at an angle a with respect to the cylinder axis. 

However, in several practical applications, in particular for the radar community, elongated objects may become of 
specific interest for these applications. Ref. [l5| provides quasi-static formulas for infinite circular dielectric cylinders 
under normal incidence, which was later extended to 2D infinite conducting cylinders, also illuminated at normal 
incidence in Ref. [23 |. These results were also preliminarily extended to oblique incidence in Ref. [30|. Moreover, in 
Ref. [3l| the results of Ref. i2J| for multi-frequency operation of spherical cloaks were extended to dielectric infinite 
cylinders. Also, recent theoretical and experimental efforts J32l |33| reported the practical realization of plasmonic 
cloaking in 2D cylindrical geometries, proposing metamaterial designs for specific polarization of interest. 

The literature on plasmonic cloaking applied to cylinders, however, has often dealt with idealized 2D geometries: 
infinite cylinders, and incident waves normal to the cylinder axis, with specific polarization properties. This assump- 
tion, common to several other cloaking techniques applied to cylinders, makes the resulting calculations quite limited 
from a practical standpoint. It may be argued, in particular, that once the angle of incidence is modified, and the 
end effects of finite-length cylinders are considered, such cloaking effects may be severely limited, if not completely 
lost. In [SJ, [33, the effects of truncation were preliminarily considered for normal incidence. 

Here we analyze all these issues in great detail, first deriving a general cloaking theory for arbitrary illumination 
of an infinite cylinder. We show that it is indeed possible to find a suitable, robust plasmonic cloaking layer in this 
scenario, which under suitable conditions may operate over a broad range of incidence angles. We corroborate these 
findings with an extensive numerical analysis of finite length and truncation effects in practical cylindrical geometries, 
studying the overall scattering reduction for different angles of incidence and several design parameters. The present 
analysis applies to idealized metamaterial cloaks with isotropic properties. We leave to future work the practical 
limitations introduced by the specific realization of engineered metamaterial cloaks. 



II. ANALYTICAL RESULTS FOR INFINITE CYLINDERS 



General Derivation 



Consider the geometry of Fig. [U depicting a circular cylinder of finite length L, radius a, permittivity e and 
permeability /x, covered by a thin conformal cylindrical cloak shell of thickness {ac~a)^ permittivity tc and permeability 
IXc- In this section, we examine the limiting case of an infinite circular cylinder (L — > 00) illuminated by an arbitrary 
plane wave at oblique angle. The general scatteri ng p roblem may be solved by expanding the impinging and scattered 
fields in cylindrical harmonics (see e.g. Refs. |36l - [38| ). which we apply to derive the scattering response as a function 
of incidence angle a. 

Without loss of generality, the incident wave may be decomposed into its transverse- magnetic (TM^) and transverse- 
electric (TE2) components, with respect to the cylinder (z) axis. By matching boundary conditions at the two radial 

Analytical results for normal 

In that special case (a = 90°), 

H, for TE. 
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interfaces, the problem may be solved exactly for each cylindrical harmonic [3 

incidence, specifically applied to the cloaking problem at hand, were derived in Ref. 

only two tangential field components are non-zero at each interface: Ez,H^ for TM polarization, and E, 

This ensures that the scattering problem is easily solved as a set of four equations, and the TE and TM Mie scattering 

coefficients are easily expressed as a function of 4 x 4 determinants, [/^^'"^^ and V^^'"^^ |32| : 
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(1) 



for cither TM or TE polarization. We have assumed here and in the following an e~"^* time dependence. The overall 
scattering cfBcicncy, defined as the ratio of the total scattering cross section normalized to the cylinder's physical 
cross section, is given by [37| : 
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As discussed in [32[, by tuning the cloak's electrical and geometrical properties, the cylinder's visibility can be 
minimized by canceling the dominant scattering orders, which is possible under the condition [/„ = 0. In the quasi- 
static limit, i.e., for very thin cylinder s, appr oximate closed-form expressions for these cloaking conditions for TE and 
TM incidence were reported in Refs. J29ll39l|. 

Scattering for oblique incidence is a more challenging problem, since polarization coupling is usually involved [38| . 
as discussed in the Appendix □ This implies that for oblique incidence purely TE or TM cylindrical waves may excite 
also TM and TE scattered waves, respectively. This feature will play an important role when analyzing the response 
of the cloak for different incidence angles. In the special circumstance for which this coupling may be absent, or 
minimized, e.g., as discussed in the Appendix, for azimuthally symmetric modes (n = 0), for conducting or high- 
contrast cylinders, or for small cross-sections, the scattering coefficients may be still expressed as in Eq. (1), where 
Un and Vn have generalized expressions: 
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where J„ and Yn are the cylindrical Bessel functions of the first and second kind of order n, J' ,Y' are their derivatives 
with respect to the argument, and i]i = \/ [i-ij ^i are the characteristic impedances of each region. The TE expressions 
may be found by electromagnetic duality: e^ ^ Uj. 

These expressions coincide with those in Ref. [S^] for normal incidence (a = 90°), as expected. The dependence 
on the angle of incidence a is encoded in the transverse wavenumbers, kf . Specifically, kf = ^fc? — ^'^, where 
/3 = kocosa is the wave number component along the cylinder axis. In particular, fej^ = fegsino;, as expected. In 
the general case, where polarization coupling will occur, Ec^. [1] should be modified to include the TM-TE coupling 
coefficients, consistent with the derivation in the Appendix and the results reported in [36]. 



B. Quasi-Static Analysis 



As in the spherical scattering scenario |15| and cylindrical scenario at normal incidence [32| , it is instructive to first 
consider the case of electrically small cylinders, i.e. kiUc <C 1. In this limit, as we discuss in the following, the analysis 
is made simpler by the fact that: (a) the dependence of Eqs. pi4p on the incidence angle is negligible; consequently, 
so is the cross-polarization effect on the cloaking conditions; (b) the scattered wave is dominated by a limited number 
of multipolar orders. The combination of these features makes the plasmonic cloaking easier to achieve, more effective 
and robust, compared to larger cylinders. This is expected, since satisfying the condition C/„ = for only few dominant 
terms in Eq. ^ is easier to achieve, and such conditions are less dependent on the angle of incidence in this regime. 
In particular, it is easy to show that for regular dielectric cylinders (e > eo, ^ = /io)j the dominant terms in Eq. ([2]) 
are c™ and c™, respectively for TM and TE excitations. For arbitrary polarization, c™ scattering dominates. This 



also applies to a conducting cylinder, whose limiting case corresponds to e — >■ —joo and /i 
magnetic cylinder (^ > fiQ, e = eg) dual considerations apply, so c™ dominates. 



in Eqs. dSU). For a 



1. Quasi-static cloaking conditions 

As discussed in the Appendix, for electrically small cylinders {kiUc <^ 1), the cross-coupling terms vanish for any 
angular order n. We may thus derive the quasi-static cloaking condition by simply taking the first-order Taylor 
expansion of the coefficients of Eq. ([T]) . This leads to simple closed- form solutions for the appropriate ratio of cloak 
to core radii to achieve cloaking in the quasi-static limit for each scattering coefficient, consistent with analogous 
formulas available in the spherical [l5| and cylindrical normal-incidence |32| scenarios: 
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For perfectly electric conducting (PEC) cylinders, a case of particular interest for cloaking applications at radio 
frequencies, we find: 
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We note that Ref. [29| misreported that there is no possible quasi-static condition for c™. The correct formula, 
shown above, results from taking the proper limits in our previous analysis for e -^ — joo and /x — > in U^^. This 
is particularly relevant for oblique incidence, since enforcing only e — ^ —joo to model a PEC boundary, as suggested 
in Ref. |29i , would not ensure zero normal component of the magnetic field on the boundary, as required in the PEC 
limit [4^. As correctly reported in [23: there is no quasi-static condition for the c™ coefficient in the PEC scenario. 

A few important points should be highlighted regarding Eqs. ^H^. First of all, Eq. ([5]) for TE polarization 
formally coincides with the derivation of Ref. [15|. The formulas here have been obtained for completely generic 
oblique incidence, generalizing the previously published result. They show that in the long wavelength limit, cloaking 
conditions are unaffected by an arbitrary variation of the angle of incidence. In this regime, the TM-TE polarization 
cross terms are of second-order (cf. Appendix ^ , thus Eq. ([5]) ensures that an electrically small magnetodielectric 
cylinder may always be cloaked using its dominant scattering order in the Taylor expansion, regardless of the angle 
of incidence. As in the spherical scenario, these quasi-static formulas split the role of permittivity and permeability 
between TM and TE polarizations. Special attention should be paid to the azimuthally symmetric modes, for which 
the role of permittivity or permeability is reversed compared to higher-order modes. Moreover, again as in the spherical 
case, the cloaking conditions depend only on the ratio Gc/a, implying that a thin homogeneous shell may be employed 
to suppress the relevant scattering orders in this quasi-static scenario. It should be stressed that Eqs. ^H^ may not 
be met by arbitrary values of the constituent parameters. Rather, they should lie in specific ranges of permittivity, 
due to the simple physical constraint on geometry ac/a > 1. 

From a practical standpoint, an electrically thin dielectric cylinder, for which the c™ coefficient is dominant, may 
always be cloaked by a thin shell with Cc < eg. Under the first of the conditions in Eq. [Sj most of the scattering 
may be suppressed, although significantly negative Cc may be required to achieve a very thin shell. By duality, an 
electrically thin magnetic cylinder will require /ic < Mo- 

For such infinite cylinders, the cloaking conditions may in general become trickier than in the spherical case, 
since the azimuthally symmetric cylindrical waves follow special cloaking conditions: they are permittivity-based for 
zeroth-order when higher order arc permeability-based, and vice versa. This is particularly relevant for the conducting 
scenario, for which the dominant c™ coefficient may not be canceled at all in the long-wavelength limit (cf. Eq. ([6])). 
It should be emphasized that it is still possible to suppress scattering in the fully dynamic scenario, since it may always 
be canceled with its dynamic expression Uq^. However, in the long- wavelength limit, the required wave number in 



the cloak considerably grows, implying that quasi-static considerations such as those used in Eq. (j6]) may not be 
applied. However, it is evident that conducting objects are significantly more challenging to cloak for cylinders than 
spheres [20] • The physical reason for this difference is evidently related to the fact that electrically thin cylinders are 
not electrically small objects. These systems are still infinite in the axial direction (for this analytical treatment), 
which implies that their scattering may not necessarily be small - above all when conduction currents may be induced 
in the z direction, as with the TM^ cylindrical wave. 

2. Validation of the quasi-static conditions: Electrically small dielectric cylinder 

To assess the validity of the previous approximate cloaking conditions in the fully dynamic scattering problem, 
we present some numerical examples of interest in this section, solved using the exact analytical theory reported 
above. Consider first the case of a dielectric cylinder with e = 3eo and normalized diameter /coQc = 0.1, covered by 
a thin uniform cloaking shell with Uc/a =1.1. Fig. [2^ shows the variation of the total scattering efRciency Cs as a 
function of the cloak permittivity and of the angle of incidence of a TM illuminating plane wave. It is immediately 
apparent that the use of negative-permittivity thin cloaks may significantly reduce the overall scattering efficiency of 
the cylinder, and in particular minimum scattering is obtained very close to the corresponding quasi-static cloaking 
condition for c™ in Eq. ([5]), which would yield Ec ~ — 8.524eo. Moreover, minimum scattering is achieved at the 
same value of negative permittivity, independent of the angle of incidence, consistent with the previous theoretical 
results, despite some polarization coupling for oblique angles. Overall scattering reduction, compared to no cloak 
(ec — eo)j is especially significant in the normal incidence case, which provides over 50 dB reduction. For oblique 
incidence, the coupling with TE coefficients affects the cloaking performance and generates residual scattering at 
the cloaking condition, but it is seen that the same (negative) value of cloak permittivity may provide significant 
scattering reduction (over 10 dB) for a wide angular range without drastic changes. 

For smaller incidence angles, we also observe excitation of a plasmonic resonance for slightly negative values of Ec- 
This is associated with the condition V^^ = 0, obtained here in the quasi-static limit for Ec = — 0.07eo, consistent with 
findings on quasi-static plasmonic resonances in layered spheres and cylinders |39j . It is evident that the resonance 
is not excited at normal incidence, but arises as soon as oblique incidence is considered, due to cross-polarization 
coupling. This explains the sharp peak for small negative values of Ec in Fig. [2^. Analogously, for low positive values 
of ec, a second small dip in the scattering efficiency appears at a = 7r/4, associated with the cloaking condition for 
c™ = in Eq. ([5]). Here, this provides the solution ec = 0.14eo. Once again, this dip does not appear at normal 
incidence, due to the lack of coupling between polarizations for a = 7r/2. To conclude this discussion, it should be 
underlined that, in absolute value, scattering is larger for normal incidence in the uncloaked scenario (ec ~ eo in the 
figure), as expected, due to the larger Ez component. On the other hand, the proper cloak design may maximally 
suppress scattering in this worst-case condition, which is an inherent property of this cloaking approach. 



Dielectric Cylinder 




-20-15-10-5 

ecAo 



10 15 20 



10 



-10 


Magnetic Cylinder 


-90° 
- 67.5° 
45° 




(b)_ 


-20 






^ -30 


' ^ 


^___^ 


r 




..----^ 


-50 


\f 


- 


-60 


- \ 


- 


-70 
on 





-25 



-20 -15 -10 

Mc/Mo 



-5 



FIG. 2: Total scattering efficiency for electrically thin infinite-length circular cylinders with fcoOc = 0.1 and Oc = 1.1a: 

(a) dielectric (e = 3eo, A* = A'o) and (b) magnetic (e = eo, ^ = 3/io); for three different angles of incidence of TM^-polarized 

radiation. In the electrically thin limit, the cloaking conditions have weak dependence on the angle of incidence. 



3. Validation of the quasi-static conditions: Electrically small magnetic cylinder 

As a second example, consider the same cylindrical geometry, but with magnetic properties {fj. = 3/io), again excited 
by a TM^ wave. All other parameters are kept the same. The cylinder's scattering efficiency is shown in Fig. [5]3. 
As anticipated above, this configuration gives much weaker scattering in the uncloaked scenario (fic = Mo)i due to its 
small size and lack of dielectric contrast with the background. In fact, the c™ coefficient is negligible in this case. 
Residual scattering is dominated by the higher-order cj^ coefficient, which may be canceled with ji^ = — 21.2/io or 
0.14/zo (note the dual behavior compared to c™ in Fig. [2fe.)- Two clear scattering dips are visible around these two 
values, which also in this scenario do not depend on the angle of incidence, consistent with Eq. [SJ The absence of 
scattering from higher-order modes and of coupling for the n = modes ensures that overall cloak performance is 
unchanged by variations of the incidence angle a, as evidenced by the different curves in Fig. [2}d. 

4- Quasi-Static Conditions: Conclusions 

It should be stressed that the above numerical analyses neglect Ohniic absorption losses in the cloak and/or core 
materials. We include these effects in the following Sections. Regardless, the plasmonic cloaking technique has 
been shown to be inherently robust against moderate absorption [20], since it is not based on a resonant effect. 
This implies that the curves in Fig. [2] would remain practically unchanged near the cloaking regions when moderate 
losses are considered. On the other hand, the large scattering peaks associated with plasmonic resonances would be 
significantly dampened by realistic losses. This and the previous considerations imply that: 

(a) The plasmonic cloaking technique may be successfully applied to infinite cylinders. 

(b) The quasi-static conditions, Eqs. (I5]|6]), hold to a very good approximation for thin cylinders, fco^ic ^ 0.1, and 
ensure dramatic scattering reduction in this limit. 

(c) Cloaking conditions for electrically thin cylinders arc very weakly dependent on a. 

Thus, cloak designs for electrically thin cylinders are robust against the angle of incidence for any scattering order n 
and polarization. We explore numerically in the following Sections how these properties are affected by considering 
thicker geometries and truncation effects. 

III. PLASMONIC CLOAKING DESIGN AND OPTIMIZATION 

For a given infinite cylinder of radius a and permittivity e, a single-layer plasmonic cloaking shell may be optimized 
by varying its two design parameters Oc (or Uc/a) and Cc, using the analytical results derived in Section |lll In general, 
it is preferable to choose thin shell thicknesses, since drastically increased cross-sections imply reduced bandwidths 
and larger sensitivities to the design parameters [231 . 

This simple cloaking design inherently limits the overall total thickness of cylinders that we may cloak, since only 
few scattering orders may be drastically suppressed independently. Magnetic properties of the cloak [l9[ as well as 
multi-layer designs [2J| may be considered to increase the available degrees of freedom and size of objects to be cloaked. 
In this work, however, for sake of simplicity we limit our interest to non-magnetic cloaks, and we consider moderate 
cross-sections of the objects to be cloaked. In particular, in the following we consider cylinders with diameters 
2a = (^, ^, ^) and relative permittivities of 3 and 10 in the dielectric scenario, in addition to a PEC cylinder. 
Since the previous theoretical discussion shows that infinite cylinders may be effectively cloaked, as anticipated, we 
consider overall lengths L of several times the diameter, which may be comparable to or larger than the wavelength 
of operation. 

Because of the TM-TE polarization coupling described in the previous Section (cf. also the Appendix J, completely 
general cloak optimization would be quite complicated, and in general depend on the incidence angle and polarization 
of excitation. We choose instead to optimize cloak response at normal incidence - which is usually the angle for 
which larger scattering is produced - then examine the response at oblique angles with numerical simulations. More 
extensive cloak optimization would involve the choice of a cost function and the application of global (e.g. genetic 
algorithms [41|) and local (e.g. Quasi- Newton) optimization techniques across all angles and polarizations, also as a 
function of the observer's position. In general, further optimization may be achieved by considering magnetic cloaks, 
as discussed above. We will analyze these aspects in the near future. 



A. Design Parameter Space: the View from Above 

To provide insight into the effects of the various parameters involved in cloak design, we examine the variation of 
scattering efficiency with ad a and tc for TE or TM plane wave normal incidence in the infinite (2D) scenario, using 
the formulation developed in the previous Section. To that end we write the scattering gain as: 
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(7) 



Cs is the scattering efficiency from Eq. ([2]), and the superscript represents the uncloaked case, calculated with 
Eqs. (HJHI), in the limit a^-^ a and e^ — > e. 

In general, the above discussions show that it is not possible to achieve large scattering reduction for both TM 
and TE illumination simultaneously with a single-layer permittivity cloak at the same frequency. For moderate cross- 
sections, however, such as those considered here, scattering is generally dominated by one of the two responses. It is 
thus best to design for maximal suppression of that polarization. Hence, we focus on TM^ waves, which dominate 
scattering from infinite dielectric and conducting cylinders in the long wavelength limit (cf. quasi-static analysis in 
the previous Section). Ideally, the cloak is designed to provide minimized scattering gain for normal incidence. 

For our calculations we used Mathematica [42| and truncated the Mie summations at order n = 5, which ensured 
convergence for the cross-sections considered here. Fig. |3] shows scattering gain contour plots for infinite cylinders of 
diameter 2a = Ao/4, relative permittivities e = 3 (Fig. [3^) and e = 10 (Fig.lSja), and TM normal incidence. The plots 
show alternating loci of resonant peaks (large Qs^ light color) and cloaking regions (near-zero Qs-, dark color). As 
expected, thinner cloaks (smaller ad a) yield larger bandwidths, reflected by wider cloaking regions. This implies that 
the ideal cloak design would utilize Ec < 0, as expected from the results in Sec. |lTl from simple physical considerations 
and analogous results in the spherical scenario [15|. The cloaking loci for large positive values of ec are associated 
with anti-resonances that arise when the shell thickness is comparable to the wavelength in the shell material. These 
anti-resonances are narrow-bandwidth and highly sensitive to the design parameters, and they necessarily lie in close 
proximity with scattering resonant peaks, which makes them not suitable for a robust cloaking design. Fig. |4] shows 
constant ad a slices of Fig. [3] to better illustrate these characteristics. These slices illustrate how ec < choices 
typically produce lower scattering gain, i.e. better cloaking, away from dangerous resonant enhancements. This is 
especially true when cloaking low-density dielectric cylinders. 

We also considered optimization of plasmonic cloaks for PEC cylinders of the same size. However, as noted in the 
quasi-static limit, here the dominant scattering contribution (via the c™ coefficient), may not be canceled in the 
long-wavelength regime. This implies that in the corresponding contour plots there would not be cloaking regions for 
thin negative-permittivity shells, as in the dielectric scenario. Moderate scattering reduction may be achieved with 
large-ec thick shells, but as outlined above this implies large sensitivity to frequency, design variations and closely- 
spaced resonant peaks, which may be excited at different incidence angles. It is thus evident that a simple permittivity 
cloak may not be sufficient to adequately cloak a PEC cylinder and, as in the spherical scenario [15|, it may require 
magnetic permeability ^c different from the background for robust scattering reduction. In the following, therefore, 
we mainly focus our design efforts on dielectric cylinders. 
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FIG. 4: The ratio Qs as a function of cloak permittivity, £2 for slices of constant cloak thickness, ac — 1.1a, in the contour 
plots of Fig. for core permittivities e = 3 (left) and e = 10 (right). 

B. Cloak Design 

The previous results show that optimal cloak configurations for dielectric cylinders are based on negative permittivity 
metamaterials and thin cloak shells. We should point out that negative values of effective permittivity may indeed 
be achieved in the microwave or THz frequency ranges using various metamaterial geometries, such as wire media 
or parallcl-platc implants |32| . and they are naturally available at larger frequencies. In particular, the parallel- 
plate implant technology may be particularly well-suited for cloaking incident TM^ waves, as it has been already 
demonstrated theoretically and verified experimentally for normal incidence at microwave frequencies j32l l33j . We 
assume here, for sake of simplicity, that the required value of effective permittivity is available for the shell geometry 
of interest and that the cloak material is isotropic. Possible anisotropy for a specific metamaterial realization, inherent 
in some proposed realizations 13^|33[, may affect cloak performance for different incidence angles, but our preliminary 
results show that for thin cloaks such effects may be minor. Thus, for simplicity we always assume idealized isotropic 
metamaterials in the following. The requirement that any passive metamaterial has a frequency dispersion ensuring 
d{ujec) / duj > is met here by assuming a Drude dispersion model of the form: 
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Here and in the numerical simulations in the next Section, we calculated the plasma frequency LOp to ensure that at 
the design frequency /g the real part of e^ yields the required value. The damping frequency 7, associated with the 
level of losses in the metamaterial, has been always assumed in the following to be 7 = 10~^LUp. This value provides 
a moderate amount of loss, comparable with practically-realizable metamaterial geometries at these frequencies. 

Table U summarizes an extensive campaign of optimizations that we performed to cloak different cylindrical ge- 
ometries, as outlined above, considering several cloaking designs. Inspecting Table I, we see that thin cylinders 
(2a = A/8), consistent with the previous quasi-static analysis, require optimized cloaks with negative permittivity in 
the case of dielectric objects, and a large positive permittivity for conducting materials. The corresponding scattering 
gain may become extremely low, since only few Mie coefficients effectively contribute to the total scattering, and 
they are properly canceled by the optimized cloak design. As a general rule of thumb, consistent with the previous 
considerations, thinner cloaks require larger values of permittivity, either positive or negative. Relatively thicker 
cloaks relax the requirements on very negative permittivity of the cloak, but, as highlighted above, larger thickness 
also tends to produce additional scattering terms, which limits overall performance. However, lower absolute values of 
negative (in the case of dielectric) or positive (for PEC) permittivity may be easier to achieve and be less sensitive to 
loss. A trade-off between cloak thickness, material reliability and overall cloaking performance should be considered. 
Larger dielectric objects may be cloaked by thicker cloaks with positive values of permittivity, since the dynamic 
nature of the wave in the cloak may produce a negative polarizability even with positive materials in larger shells. 
However, the overall scattering reduction is less dramatic; residual scattering is indeed expected for larger objects. 
In particular, as discussed above, conducting cylinders are most challenging to be cloaked against TM waves, due to 
the axial conduction currents induced on their surface. Scattering reduction of around 3 dB, however, may still be 
achieved even for conducting cylinders. Moreover, large values of core permittivity e arc harder to cloak, and the 
residual scattering is larger than for a lower-permittivity core of the same size. 

Our cloak designs in Table U were optimized for TM normal incidence, the scenario with maximum scattering for 
the considered objects. At oblique incidence, however, TM-TE coupling excites additional scattering modes, which 
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TABLE I: Parameters from the cloak optimization procedure described in Sec. lIIIBl for different objects and cloak thicknesses. 
Qa corresponds to the scattering gain for normal-incidence TM2 illumination. 



may affect the overall response, in particular for thicker objects and moving toward grazing incidence. We discuss 
these features in detail in the next Section. 

Fig. [5] considers the typical frequency dependence of cloak performance on wave polarization. The left panel shows 
the uncloaked scattering widths for an infinite dielectric cylinder with 2a ~ Ao/4 and e = 3 for TM and TE normally- 
incident plane waves. As expected, the TM polarization is more strongly scattered, thus it is of more relevance in 
designing an optimized cloak. However, despite the much smaller response to TE waves, it is interesting to analyze the 
behavior of a TM-optimized cloak to TE polarization, since a portion of the TM scattered power may be coupled to 
TE waves for TM oblique incidence, in particular for smaller incidence angles. The right panel shows the normalized 
scattering gain for TM and TE incidence, considering a conformal shell cloak with Oc = 1.1a and ec — —13.55 at design 
frequency /o with Drude dispersion, which coincides with the optimized design at /o for TM normal incidence. We 
see how the cloak is very effective over a relatively broad range of frequencies around /o for TM incidence, despite the 
natural dispersion of the cloak material. However, this same cloak does not provide scattering reduction at the same 
frequency for impinging TE waves, supporting instead a small cloaking dip at lower frequencies. Although this per se 
is not a major issue, since the absolute value of the TE scattering is much smaller (left panel), and the quantities in 
the right panel are normalized to the corresponding uncloaked scattering, this issue may cause an overall dependence 
of the cloaking response to the incidence angle for TM excitation, due to the TM-TE polarization coupling. We verify 
these issues in the next Section, which considers the performance of finite-length cylinders illuminated at oblique 
incidence. 
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FIG. 5: Uncloaked scattering efficiency (a) and expected scattering gain (b), Ca/C^, as a function of relative frequency, for an 
infinite cylinder with 2a = Ao/4, e = 3 and for Oc = 1.1a (see Table|l|. Optimized for TM waves at normal incidence, the cloak 
will perform differently for TE illumination. 
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IV. NUMERICAL RESULTS FOR FINITE-LENGTH CYLINDERS AND OBLIQUE INCIDENCE 



After having considered the optimization of cloaks for dielectric and conducting infinite cylinders, we consider in 
this Section the effects of finite length and varying angle of incidence on the performance of cylindrical plasmonic 
cloaks, using full- wave numerical simulations [43|. To keep the design simple, in our numerical simulations we simply 
truncated the cloaked cylinder, as optimized in the previous section, consistent with the geometry in Fig. [1] leaving 
the ends uncovered. This is expected to cause additional scattering, especially for more oblique incidence (smaller 
a). However, since this cloaking technique is based on an integral effect [20|, for moderate cross-sections these effects 
are not expected to significantly deteriorate the overall cloaking performance. This intuition is confirmed by the 
following numerical results. Improved performance may be obtained by locally tailoring the cloak around the edges, 
and covering the ends, as described for conical geometries in [28|. The small improvements achievable with this 
fine-tuning of the cloak design are in any case not relevant to the following general discussion. 

A. Dielectric Cylinder at Normal Incidence 

Consider first the dielectric cylinder of Fig. [2^ (e = 3eo, fJ. = fJ.o, diameter 2a — 2.5 cm), but with finite length 
i = 10 cm. Our objective is to reduce its scattering cross-section at the design frequency /o = 3 GHz, for which 
this object has cross-sectional width 2a = Ao/4 and length L = Aq; Aq = 10 cm is the free-space wavelength at 
frequency /q. In this frequency range, the object may produce significant scattering, since its length and electrical 
size are comparable to the incident wavelength. Scattering at normal incidence is dominated by the TM^ polarization 
(cf. Sec. ini, so we excite the cylinder with an impinging TM wave, and use a cloak with permittivity e = — 13.55eo 
{ujp = 11.44 GHz) at /q and thickness Uc = 1.1a, obtained in Sec. IIIII Note the relatively large negative permittivity 
required by such thin layer. For a thicker cloak, jed may be reduced, but at the price of possible excitation of 
higher-order cylindrical harmonics, which may slightly deteriorate the cloak performance. 

Fig. [6^ plots the total 3D radar cross section (RCS) l3g, obtained by integrating the scattering cross section over 
all angles, calculated using commercial software based on the finite-integration technique [43f. It compares the 
cloaked (solid green curve) and uncloaked (red dashed curve) scenarios for normal incidence {a = 7r/2). Despite 
the cylinder's finite electrical length and its non-negligible cross-section, scattering at the design frequency may 
be significantly suppressed with a proper cloak design. As expected, the scattering is a minimum at /o, and the 
bandwidth of operation is moderately large, as is common for plasmonic cloaks [20[ . Significant scattering reduction 
is achieved over a fractional bandwidth of over 30%. Fig. |B}d plots the corresponding scattering gain, defined as 
the ratio between the total RCS of the cloaked and uncloaked geometries. The curves effectively show the overall 
scattering reduction achieved by adding the cloak to the bare cylinder for different frequencies. The plot compares the 
full- wave simulations for finite-length cylinders calculated using [43|, with the infinite-cylinder analytical formulation 
derived in the previous section. The overall scattering reduction is « 15 dB at /o, and the analytical result for an 
infinite cylinder agrees extremely well with the full-wave simulation for a finite cylinder, despite the truncation effects. 
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FIG. 6: (a) Frequency dispersion of the total RCS for a dielectric cylinder with e = 3eo, n — ^o, 2a — 2.5 cm and finite length 
L = 10 cm, comparing uncloaked v. cloaked (as optimized in Section llll|) with ac = 1.1a and Re[ec]= — 13.55eo using a Drude 
model at /o = 3 GHz for TM illumination at normal incidence {ujp = 11.44 GHz), (b) The corresponding scattering gain 
(ratio of cloaked to uncloaked RCS), compared with analytical results obtained for an infinite cylinder with otherwise identical 
geometry. The full- wave simulations were performed with [431 • 
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TABLE II: Scattering gain at frequency /o for the cloak designs of Table|T]with 10% thickness (qc = 1.1a), and TM^ normal 
incidence. We compare the theoretical results for infinite cylinders (Qa,th) with the numerical full-wave simulations for truncated 
cylinders with L = \o {Qs.si,n)- 

At lower frequencies, < 2 GHz, the plasmonic nature of the cloak and its negative permittivity induces a plasmonic 
wave along the cylinder, which actually increases the overall scattering compared to the uncloaked scenario. This is 
predicted in both the infinite and finite-length cases, and readily observed in Fig. |B}d. There, the contribution is even 
stronger for the finite-length truncated cylinder, due to plasmonic radiation at the ends. 

Table HIl summarizes the calculated overall suppression at /o for the dielectric-core Uc = 1.1a cloak designs of Table |I] 
illuminated by TM waves at normal incidence. In particular, the table compares the simulated scattering gain for 
finite lengths Cs,sym with the theoretical value predicted analytically for infinite cylinders Cs,th- We observe excellent 
agreement between predicted results from the analytical formulas for the infinite cylinder, and full-wave numerical 
simulations for the truncated geometries. This implies that truncation effects are not significant in several realistic 
setups when cloaking dielectric cylinders with the plasmonic cloaking technique. 

B. Varying Angle of Incidence 

Fig. [7] shows the variation of scattering gain for the cylinder of Fig. [5] varying the incidence angle a, as depicted 
in Fig. [TJ Even for small (near grazing) incidence angles, the cloaking effects are barely different from the normal- 
incidence case. For smaller values of a, cloaking is somewhat reduced at /g, due to the excitation of TE scattered 
waves via cross-polarization coupling, as discussed in Sec.|lTl These results demonstrate strong agreement between the 
analytical calculations for infinite cylinders and the numerical simulations for finite L, although some expected minor 
deviation does appear at small angles, due to end effects. We emphasized that, although the normalized scattering 
gain is reduced less for smaller angles, the overall RCS is significantly smaller in these cases, since the electric field 
component along the cylinder is comparatively shorter, consistent with the discussion in the previous Section. The 
results of Fig. [7] demonstrate convincingly that plasmonic cloaking may be applied to finite- length dielectric cylinders 
at oblique excitation, although at small a one may be required to cloak the TE contribution separately for improved 
performance. This is consistent with our findings in Fig. [5] 
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FIG. 7: Scattering gain for the cylinder of Fig. [51 for various angles of incidence a. (a) Results from numerical simulation, (b) 
Numerical simulation (solid) v. analytical infinite-cylinder results (dashed) for select angles. 
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FIG. 8: Left: Comparison of scattering gain for two diameters of dielectric cylinders, as labeled, illuminated at normal incidence 
with a TM polarized wave. Right: Comparison of scattering gain under normal-incidence TM illumination for a cylinder with 
same geometry as in Fig. |B] but with core permittivities e/eo = 3, 10 as labeled. The calculations used the optimized cloak 
designs as described in Table IT] 

C. Varying Diameter 

The left panel of Fig. [5] shows the scattering gain under normal-incidence TM^ illumination for the cylinder of 
Sec. IIV Al (cf. Fig. [6]) next to that for a cylinder with twice that diameter, 2a = Ao/2 (cf. Table U for cloak 
parameters). The figure also compares these results with the infinite-cylinder analytical result of Sec.|lTl End effects 
are negligible except at low /, and the analytical results match very well the numerical simulations - especially at 
frequencies where plasmonic resonances are not excited. Cloaking is effective over a relatively broad bandwidth, 
even for the thicker cylinder. This implies that the designed cloaks may be effective over a relatively broad range of 
object size. The excellent agreement between the analytical curves for infinite cylinders and the simulation results for 
truncated geometries indicate that truncation effects are negligible for cloaking performance and suggest the use of 
the previously derived analytical formulas for fast cloak optimization. A simple permittivity cloak, as considered here, 
is effective for significant scattering reduction for cylinders with diameters of the order of the wavelength. Thicker 
cylinders may require the use of multilayer and/or magnetic cloaks, as discussed in Ref. [19[ for spherical geometries. 



D. Varying Object Permittivity 



The right panel of Fig. [5] shows the scattering gain for the cylinder of Sec. IIV Al (cf. Fig. IH]), for the same dielectric 
and also with a denser core, e = lOeo (see Table |T|. In the large-e limit, this comparison underscores some of the 
challenges that may be involved in cloaking a conducting object, which would coincide with a dielectric core in the 
limit of very large permittivity. There is significant scattering reduction around the design frequency /o also in the 
denser scenario, but this example simultaneously supports even stronger cloaking at a lower frequencies, / ^ 2 GHz. 
This effect, predicted by the analytical results for infinite cylinders, is associated with the frequency dispersion of 
the cloak, which matches the cloaking condition for the coefficient at lower frequencies (for more negative values). 
Effectively, this scenario presents a coincidental cloaking effect at another frequency. This does not imply that one 
could tune the second suppression arbitrarily, since its position depends on the natural metamaterial dispersion. 
However, one could tune its separation from /o to a limited degree if the cloak thickness were not fixed at Oc ~ 1.1a. 



E. Varying Cylinder Length 



Fig. |9] analyzes in more detail the effects of truncation for the case of Fig. |6] in Sec. IIV Al by considering different 
truncation lengths for several incidence angles: a ~ i^ JI^'k j'i^'K j& for the three panels, respectively. Each compares 
different lengths, including the analytical result for infinite cylinders. The cloaking effect is indeed robust against 
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FIG. 9: Scattering gain for the cylinder of Fig. [6] and variations in length by a factor of two, for three different ang 
incidence under TM excitation. The cloaking design is robust even for very short lengths, almost up to L ~ 2a. 
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variations in incidence angle, and truncation effects are quite moderate, as the scattering gain follows the line calculated 
analytically in the infinite cylinder geometry. Even for the shortest length (L = 5 cm, a 2:1 aspect ratio), the full-wave 
simulation still follows the infinite- length case with surprising accuracy, above all near /q. In this case, scattering is 
characterized by small resonant peaks at lower frequencies, associated with longitudinal resonances due to the finite 
length. The results are nevertheless extremely encouraging, in particular for shorter cylinders, for which the cloaking 
design works extremely well in the frequency range of interest. 

F. Near- and Far-Field Plots 

Figs. [TOlfTT] illustrate the full- wave near- and far-field numerical results for the cylinder of Sec. IIV Al (Fig. [6|), 
illuminated by a TM wave impinging at an angle a = tt/3. The various figures compare in panel (b) the cloaked 
configuration at the design frequency /o with the uncloaked case (a). Consistent with Fig. [5^, the overall calculated 
RCS reduction at the central frequency is -14.2 dB, and the following figures depict the effective functionality of the 
cloak in its near- and far-field regions. We have chosen a = 7r/3 to indicate the performance of the cloak to oblique 
excitation; similar considerations apply for any other angle of incidence, which we have verified via detailed study. 

Fig. [TOl shows the (normal) magnetic field distribution on the E plane (snapshot in time), comparing the uncloaked 
scenario (a) with the cloaked one (b). Several interesting features are noteworthy. In the uncloaked case, the wave 
penetrates the dielectric rod and experiences a wavelength shortening that effectively distorts the planar wave fronts on 
the back of the cylinder, producing significant shadow and scattered fields all around the object. The thin plasmonic 
metamaterial shell is able to re-establish the proper planar fronts just outside the cloak (b), ensuring reduced scattering 
and suppressed visibility for an outside observer positioned anywhere in the near- or far-field of the object. It should be 
noted that this effect is obtained for oblique incidence, and the cylinder ends are uncloaked. Additional improvement 
may be achieved by proper cloaking of these terminations. We may also easily observe how the plasmonic layer 
supports surface-plasmon waves traveling along the shell, as expected due to its negative permittivity. It is interesting 
to observe how these waves effectively cancel the residual scattering and restore the phase fronts to almost exactly 
match those of the original plane wave, had it traveled through free space instead. 

Fig. [TT] illustrates the far-field radiation patterns for this same cylinder, comparing on the same scale scattering from 
the uncloaked (a) and cloaked (b) objects, showing drastic suppression of the bistatic RCS at all angles. Uncloaked 
scattering exhibits itself mainly as a shadow on the cylinder backside; various higher-order scattering harmonics 
contribute to this residual scattering pattern. Most of the scattering is suppressed by the cloak. Panel (c) shows the 
cloaked residual scattering on a much smaller scale: as expected, scattering is not identically zero, and small lobes, 
associated with higher-order (and more directive) cylindrical harmonics not completely suppressed, are still present. 
Their relevance, however, is very limited compared to the original scattering levels. 

V. CONCLUSIONS 

We have presented an extensive investigation of the application of the plasmonic cloaking technique to circular 
cylinders illuminated by plane waves of arbitrary polarization and angle of incidence. We have derived analytical 
formulas for the general oblique-angle scenario, and showed that in the electrically-thin limit there is no angular 
dependence on the cloaking response, i.e. the design formulas are independent on the angle of incidence. 

To study the characteristics of cylindrical cloaks, we have designed a set of TM^-optimized cloaks for a wide 
range of parameters and cylinders of interest, using the normal-incidence analytical formulas with realistic losses and 
metamaterial frequency dispersion implemented in a Drude model. For the cloak design, we have focused on TM 
polarization because it dominates the scattering of moderately thick dielectric and conducting cylinders, which is 
of interest for several applications within the radar community. The optimized cloaks' drastic scattering reduction 
was corroborated with full-wave numerical simulation, taking into account variations of the angle of incidence, core 
permittivity, cylinder diameter, and most importantly truncation effects due to finite length. 

We have found that, as predicted by the analytical formulas presented here, for elongated objects with diameters 
up to one- half the wavelength (for which our cloaking technique is most effective), but with length comparable with 
the wavelength of operation, a simple one-layer permittivity cloak is very effective, providing significant scattering 
reduction highly robust to variations in the angle of incidence. Performance is slightly weakened at near-grazing angles, 
for which TM-TE polarization coupling partially affects the overall performance of a single-layer cloak optimized for 
TM polarization alone. We note, however, that scattering approaching grazing angles is the weakest in absolute value, 
thus less important for achieving overall scattering reduction. 

The analytical theory developed here for infinite cylinders has been strongly corroborated by our numerical simu- 
lations for finite lengths, implying that the truncation effects do not significantly perturb the cloaking effect. We are 
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FIG. 10: Contour plots (snapshot in time) of the magnetic field on the E plane of polarization for a TMz wave at a = n/3 
(origin top right corner) for the object of Fig. [6] (a) Uncloaked; (b) the Oc = 1.1a cloak. Severe uncloaked wavefront distortions 
are almost totally restored by the thin cloak. Cloak interface nodes (right) are associated with plasmonic surface waves. 

currently working on the extension of these concepts to multi-layered metamaterial cloaks for suppression of multiple 
scattering coefficients, and contemporary suppression of TE and TM scattered waves, which may increase the size of 
the cloaked objects and the angular range of operation. We are also currently pursuing an experimental realization 
of the plasmonic cloaking concept at radio frequencies for practical finite cylinders and cloaks. 



(a) uncloaked cylinder 



(b) cloaked cylinder 



(c) cloaked cylinder (zoomed in) 




FIG. 11: Far-field scattering patterns in the (a) uncloaked and (b) cloaked (oc = 1.1a) cylinders of Fig. [6] plotted on the same 
scale; (c) is an enlargement of (b), to show the residual detail of the scattering pattern, dominated by higher-order scattering 
modes. Panels (a) and (b) demonstrate the dramatic scattering reduction from a properly tuned plasmonic cloak. 
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Appendix: Analytical Solution of the Scattering from a Two-Layer Circular Cylinder for Oblique Incidence 

Consider the geometry of Fig. [Tl consisting of a circular cylinder of radius a, permittivity e and permeability /i, 
covered by a thin conformal cylindrical cloak shell of thickness (oc — a), permittivity e^ and permeability ^c- When 
excited by an impinging infinite plane wave of arbitrary incidence angle and polarization, the general scattering 
problem may be solved in the limit of long cylinders (L — > oo) by expanding the impinging and scattered fields 
in terms of cylindrical harmonics J36l438{ . The problem may be split into the orthogonal polarizations with _E-field 
transverse to the cylinder axis (TE^) and iJ-field transverse (TMz). 

For a TMz plane wave impinging at an angle a from the cylinder axis z, the incoming magnetic field may be written 
as: 

H, = ^e-'^^e-'''»^x , (A.l) 

770 

where Eq is the electric field amplitude, 770 and fcg are the free-space characteristic impedance and wave number, 
/3 = kocosa is the wave number component along the cylinder axis, kg = ^/k^ — /3^ = fcosina is the transverse 
component of the wave number, and the corresponding impinging electric field may be calculated using the curl 
Maxwell's equations. Expanding in cylindrical waves, we may write electric and magnetic fields as: 

^ 1 d^u™ . 1 9^u™ - , TM^ 



where 



u™ 



Eo t-'\siiia) Jn{klp) e*"-^ e""^^ (A.3) 



and J„ is the cylindrical Besscl function of order n. Analogous expressions for TE^ waves may be derived using 
duality. 

Using the orthogonality of cylindrical waves, the boundary conditions at the radial interfaces may be met by 
assuming the existence of transmitted cylindrical waves in the two dielectric regions and a scattered wave in free- 
space, which may be written consistently with Eq. (jA.2[) . using the scalar potentials: 

u™ = r^'Eo e"-^ e-'^" sina cJ^Mkfp) p<a, 

u™ = z-^Eo e^"^ e-'^^ sina [c™ J„(fcjp) + c™Y„{k:[p)] a < p < a, , (A.4) 



,TM 



I tL^e ^ e '^ sma c^, „ ii„ \Kq p) p > ac ■ 



for the fields induced in the core region, shell, and for the the scattered field, respectively. The core and cloak regions 
are labeled "1" and "2" , respectively, and "s" repres ents the scattered wave outside the cloak; ki are the relevant 
transverse wave numbers for each region and kj = yjk^ — (3"^ . Jn and ¥„ are the cylindrical Bessel functions of the 

first and second kind, for incoming and outgoing waves, while if„ is the cylindrical Hankel function. The complex 
scattering coefficients are the unknowns. Analogous equations may be written for TE^ waves. 

For normal incidence (a = 7r/2), for PEC objects, or for the azimuthally symmetric mode n = 0, the problem is 
easily solved by matching the two non-zero tangential field components at each radial interface: 

Es,z{4>) = ^2,2(0) P = Ic, 

E2,z(4>) = ^\,z{4,) P = a, 

Hs,4>{z) = ^2,,/)(z) P = Ic, 

H2M^) == ^iM^) P = a> (A-5) 

for TMz(TE2) polarization. This yields the familiar rank-four determinant expressions reported in Eqs. pi4p . Consis- 
tent expressions for the TEz coefficients may be obtained by applying duality. In the general case of oblique incidence 
on a dielectric cylinder, however, the higher-order n > modes are characterized by all four independent tangential 
components of the fields at each interface, which may not be matched independently for each TE^ or TM^ harmonic. 
The boundary conditions may instead be met, as derived in [38[, by linearly combining TE^ and TM^ harmonics of 
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same order n > 0. The corresponding solution of an eight-by-eight system of equations, derived in [S^, provides the 
exact expression of the scattering coefficients in the general case, to be used in Eq. ( [5]) to derive the total scattering 
width of the cylinder. This form of polarization coupling is inherently associated with the asymmetry introduced by 
the oblique excitation for higher-order modes, and it may be avoided only in the case of conducting objects or for 
n = 0. In the quasi-static limit, however, the dependence of the arguments of the Bessel functions in Eas. lA.4l on the 
transverse wave number kc is negligible as well; therefore this form of cross-polarization coupling is negligible. In this 
limit, Eqs. (|3]|31) may be used for any angle of incidence, at the basis of the derivation of Eqs. ([SUB]), which do not 
depend on a. 
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